Polarization angular spectra of Galactic synchrotron emission on 

arcminute scales 
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ABSTRACT 

We study the angular power spectra of the polarized component of the Galactic synchrotron 
emission in the 28-deg 2 Test Region of the Southern Galactic Plane Survey at 1.4 GHz. These 
data were obtained by the Australia Telescope Compact Array and allow us to investigate angular 
power spectra down to arcminute scales. We find that, at this frequency, the polarization spectra 
for E- and .B-modes seem to be affected by Faraday rotation produced in compact foreground 
screens. A different behavior is shown by the angular spectrum of the polarized intensity PI = 
y/Q 2 + U 2 . This is well fitted by a power law (Cpn oc t apl ) with slope ~ 1.7, which agrees with 
higher frequency results and can probably be more confidently extrapolated to the cosmological 
window. 

Subject headings: Radio continuum: ISM; Cosmic microwave background; Polarization; Methods: sta- 
tistical 



1. Introduction and main results 

In recent years the measurement of the Cos- 
mic Microwave Background (CMB) polarization 
has become one of the major aims of a large num- 
ber of planned experiments. Its detection is how- 
ever a technological challenge: so far, we have only 
upper limits on the CMB polarization level [see 
Staggs et al. (1999) for a review and the recent 
measurements by PIQUE (Hedman et al. 2001) 
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and POLAR (Keating et al. 2001) experiments]. 
Since the CMB polarization signal is expected to 
be less than 10% of the temperature anisotropies, 
instrumental sensitivities of a few /iK or less are 
required. These will be probably reached by the 
forthcoming experiments: three space missions, 
MAP (Wright 1999), Planck (De Zotti et al. 1999) 
and SPOrt (which is completely devoted to the 
study of sky polarized emissions; Carretti et al. 
2002) will measure the polarization on nearly the 
full sky, while several ground-based or balloon- 
borne experiments are planned to observe small 
sky areas with high spatial resolution [for instance, 
AMIBA (Kesteven et al. 2002), BOOMERanG 
2K2 (Masi et al. 2002) and BaR-SPOrt (Zannoni 
et al. 2002); see De Zotti (2002) for a short re- 
view] . 

The possibility of extracting information on 
cosmological parameters from CMB experiments 
is strictly related to the computation of the angu- 
lar power spectra (APS), which in the presence of 
Gaussian statistics give a complete statistical de- 
scription of the CMB emission. Although for the 
temperature fluctuations the APS is easily defined 
through the ordinary spherical-harmonic expan- 
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sion, for polarization we need two different com- 
ponents, the "electric" (E) and "magnetic" (B) 
modes, in order to describe the polarization inten- 
sity and orientation (Kamionkowski, Kosowsky & 
Stebbins 1997, Zaldarriagafe Seljak 1997). In sec- 
tion 2.2 we will discuss the definition of the polar- 
ization APS, with particular attention to the small 
scale limit. 

The detection of the CMB polarization is con- 
strained by the presence of foreground emissions. 
Different techniques have been worked out to sepa- 
rate the cosmological signal from the Galactic and 
extra-galactic emissions; all the methods exploit 
the differences in the frequency and spatial be- 
haviors. For this reason, the analysis of the APS 
has become a common tool for studying the fore- 
grounds: in fact, its knowledge allows us to esti- 
mate the foreground contamination to the CMB 
signal at different frequencies and angular scales 
(or equivalently, spherical-harmonic index £) . For 
the total intensity emission, information on the 
spatial properties of foregrounds is limited only to 
a small interval of frequencies and angular scales 
(see, e.g., Tegmark et al. 2000), and it is com- 
pletely unsatisfactory between 20 and 90 GHz. 
For the polarization, the situation is even worse 
because of the lack of high resolution surveys cov- 
ering large areas. 

In this paper we consider the synchrotron emis- 
sion, which is intrinsically highly polarized and is 
expected to be the dominant foreground at low 
frequencies. Several authors have estimated the 
APS of the polarized synchrotron using data ob- 
tained with different resolutions and from limited 
sky regions at various latitudes (for a review of 
the up-to-date surveys, see Tucci ct al. 2000). 
On small sky patches, Fourier analysis has been 
applied to the Stokes parameters Q and U (from 
which APS of the E-, S-modes can be computed) 
and to the polarized intensity PI = \/Q 2 + U 2 
(Tucci et al. 2000; Tucci et al. 2001; Bruscoli et 
al. 2002). Scalar and spin- weighted harmonic ex- 
pansions have been respectively used for PI (Bac- 
cigalupi et al. 2001; Giardino et al. 2002) and 
for Q and U (Giardino ct al. 2002). These com- 
putations are not equivalent in an important re- 
spect. Fourier analysis, being adequate for small 
sky patches, necessarily provides local effective 
spectra which for highly non-Gaussian fields (such 
as the Galactic synchrotron distribution) may be 



widely different from the global angular spectra. 
In practice such local spectra may be more im- 
portant than the global ones for the separation 
of Galactic foregrounds from CMB. Further, wc 
note that the sum Cpe = Cei + Cse, which is the 
quantity usually considered in CMB analyses, and 
Cpie provide different information. This impor- 
tant point has been noticed in Tucci et al. (2002), 
and will be discussed in Section 2.2. 

Here we observe that it helps to explain some 
discrepancies appearing in the literature. Esti- 
mates of the synchrotron polarization spectra were 
performed from the Parkes surveys of the South- 
ern (Duncan et al. 1995, 1997, hereafter D97) 
and Northern (Duncan et al. 1999, hereafter D99) 
Galactic Plane. For both surveys, sampling more 
than half of the Galactic Equator, the APS are 
nearly independent of longitude, and could be 
modeled by power laws with slopes o>e — <xb = 
1.4 -7- 1.5 and apj = 1.6 -j- 1.8 in the £-rangc 
100 ~7" 800 [Tucci ct al. 2000, 2001, Baccigalupi et 
al. 2001, Bruscoli ct al. 2002, Giardino et al. 2002 
(although the latter obtained apj = 2.37 ± 0.21 
in the Grange 40 250 from D97 data)]. Out 
of the Galactic Plane five patches are available 
at intermediate latitudes, 5° < \b\ < 20°, from 
the survey by Uyaniker et al. (1999). The APS 
vary significantly there, with slopes ranging from 
1 to ~ 2.5 (Baccigalupi et al. 2001, Bruscoli et 
al. 2002). At latitudes far from the Galactic plane 
the only available information comes from the low- 
resolution survey of Brouw & Spoelstra (1976), 
covering about 40% of the sky at five frequencies 
in the range 408 -j- 1411 MHz. From this survey, 
Bruscoli et al. (2002) found values of cxe,b, pi be- 
tween 1 -r- 2 at scales £ < 100, in agreement with 
the results from higher resolution data. Using the 
same survey, Baccigalupi et al. (2001) studied the 
PI field and found steeper spectra, with api ~ 3. 

The present paper extends the analysis of the 
APS of the polarized Galactic synchrotron to ar- 
cminute scales (£ = 10 3 -j- 10 4 ), i.e., to angular 
scales smaller than in previous works by nearly 
one order of magnitude. The study of the syn- 
chrotron contribution on these scales is relevant 
for CMB observations. In fact, the angular scales 
300 £ < 2000 are expected to be those where 
CMB should exhibit the highest level of polarized 
signal. Moreover, at £ > 3000 non-linear effects 
on CMB become important, producing polarized 
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signal stronger than the primary spectrum. These 
include the Vishniac, patchy reionization and ki- 
netic Sunyaev-Zeldovich effects (Hu 2000, Liu et 
al. 2001).' 

We make use of high-resolution polarization 
data taken from a test region for the Southern 
Galactic Plane Survey (McClure-Griffiths et al. 
2001, Gaensler et al. 2001, hereafter G01) con- 
sisting of 1.4-GHz observations carried out with 
the Australia Telescope Compact Array (ATCA). 
We find that the E and B spectra can be well ap- 
proximated by a power law at 600 ^ I & 6000, 
with a steep slope ~ 2.7 2.9. Moreover, 

we compute the spectrum of the polarized inten- 
sity PI and find that it is remarkably different 
from the above spectra, following a power law with 
a pi ~ 1.7 on the whole ^-range. We compare the 
ATCA APS with the APS computed in the corre- 
sponding patch from the 2.4 GHz Parkes survey. 
We find a noticeable agreement for Cp/^, but not 
for Ce,b,i- The different slopes found in the APS 
are interpreted in section 2.4 as due to effects of 
Faraday rotation produced by foreground screens. 
No evidence is found for a contribution of extra- 
galactic point sources. 

The behaviors of synchrotron APS on ar- 
cminute scales at GHz frequencies may be in- 
teresting for information on Galactic structure 
that spectra contain. They can tell us, in fact, 
about both the magnctohydrodynamic turbulence 
in the emitting region and the electron density 
fluctuations in the intervening medium (shown by 
Faraday rotation). In this connection, the com- 
parison between 1.4 GHz and 2.4 GHz data may 
be useful to separate the transverse structure of 
the magnetic field in the emitting region and the 
longitudinal field in the foreground screens. This 
point is open to future studies. 

2. Data analysis 

2.1. The Test Region of the Southern 
Galactic Plane Survey (SGPS) 

The data we use were obtained with the Aus- 
tralia Telescope Compact Array (ATCA), an inter- 
ferometer located near Narrabi, NSW, Australia 
and consisting of five movable 22-m antennas on 
a 3-km track (a sixth fixed antenna was not used 
in these observations). The data presented here 
were derived from observations in multiple array 



configurations (see McClure-Griffiths et al. 2001), 
resulting in virtually complete coverage visibilities 
between baselines of 31 m and 765 m. Obser- 
vations were made in the ATCA's multi-channel 
continuum mode, resulting in 9 spectral channels 
spread across a total bandwidth of 96 MHz, with 
a center frequency of 1384 MHz (see G01 for de- 
tails). There is significant Faraday rotation across 
the 96-MHz bandwidth of these data, which would 
considerably complicate our analysis. Therefore in 
the work presented here, we have only utilized im- 
ages from a single spectral channel corresponding 
to a center frequency of 1404 MHz. Faraday ro- 
tation across the 8-MHz bandwidth of this chan- 
nel is negligible. The ATCA observations covered 
a small portion of the sky, with Galactic coordi- 
nates 325?5 < I < 332?5, -0?5 < b < 3?5. This 
area was a pilot survey for the recently completed 
Southern Galactic Plane Survey (SGPS), the anal- 
ysis of which is in progress. 

Maps of the Stokes parameters Q and U were 
derived from the interferometric visibility data us- 
ing the techniques described by G01. Q and U 
images were smoothed with a gaussian beam of 
FWHM 87" x 67" and reprojected into Galactic 
coordinates. The sensitivity of the images was 
~ 1.5 mJy beam -1 , except in a strip of width 0?5 
around the edges of the field for which the sensi- 
tivity was ~ 3mJy beam -1 . 

Interferometric observations are sensitive only 
to a limited range of spatial scales. The largest 
scale that can be detected is usually determined 
by the shortest antenna spacing. In the ATCA 
case this spacing is 31 m, which corresponds to 
an angular scale of ~ 25'. Although the mosaic- 
ing process employed here allows information on 
larger scales to be recovered (Ekers & Rots 1979; 
Cornwell 1988), it does not do so with uniform sen- 
sitivity up to the largest scale sampled. Therefore, 
we here only consider angular power spectra to be 
reliable at I ^ 600, corresponding to an angular 
scale of ~ 20 arcmin. 

The distribution of linearly polarized intensity 
from the Test Region is shown in Fig. 1. The most 
prominent structure is a bright region which ex- 
tends across the longitude range 327° < I 331°, 
at latitude ~ +1°. A "spur" extends away from 
this region at I ~ 328 ?5. At the edges of the 
survey two areas with very low polarized emis- 
sion (the "voids" discussed by G01) are notice- 
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able: the first of these is approximately centered 
in (331 ?8, 1?2), the second one between 327° <, 
I < 329 ?5 at latitude - 0°. 

A wealth of small-scale structures are found in 
the diffuse polarized emission. As discussed in 
detail by G01, some of these structures are pro- 
duced by a non-uniform distribution of magnetic 
fields in the emitting regions, while in other cases 
small-scale structures in Q and U have been in- 
duced along the line of sight as a result of Fara- 
day rotation in compact clumps of ionized gas. 
Many sets of polarimetric observations have now 
identified this effect, and have attributed it to in- 
homogeneities in the interstellar medium in the 
form of clouds or filaments of sizes from 1 pc to 
100 pc (Wieringa et al. 1993; Gray et al. 1999; 
Haverkorn, Katgert & de Bruyn 2000). 

It is interesting to compare the image in Fig. 
1 with the same area of the sky observed by the 
Parkes telescope at 2.4 GHz and at the resolution 
of 10 '.4 (see Fig. 6 in G01). We expect that polar- 
ized structures induced by Faraday rotation to be 
poorly correlated between 1.4 GHz and 2.4 GHz, 
because of the strong dependence of Faraday ro- 
tation and depolarization on both frequency and 
angular resolution. On the other hand, we expect 
a good correspondence between the two data-sets 
for polarized structures which are intrinsic to the 
emitting sources. This is probably the case for the 
bright and extended emission at 327° < I ^ 331°. 
Areas of fainter polarized emission seen at 1.4 
GHz instead seem to be poorly correlated with the 
2.4-GHz data, indicating the effects of foreground 
Faraday rotation in those regions. 

2.2. The polarization angular power spec- 
tra 

The CMB temperature anisotropics are usu- 
ally expanded in spherical harmonics, AT/T(n) = 
^Zi m aimYtm{n). Assuming Gaussian statistics, 
all information in the CMB is contained in the 
angular power spectrum, defined as Ct — {\ae m \ 2 ) 7 
where ( ) denotes an average over all possible re- 
alizations of the sky. The best estimator for the 
APS is 

^ = ^ttEk™i 2 - (1) 

m 

The definition of the APS for the CMB polar- 
ization is more complex. A linearly polarized 



wave is described by the two Stokes parameters 
Q and U, which give us the polarized intensity, 
PI = \J Q 2 + U 2 , and the polarization direction 
angle <p = 0.5 atan([//Q). The value of <f> depends 
obviously on the choice of the coordinate system; 
in particular, the Q and U parameters change if 
we rotate the coordinate system. 

The standard approach to define rotationally 
invariant APS from Q and U is provided by Zal- 
darriaga & Seljak (1997) (see also Kamionkowski, 
Kosowsky & Stebbins 1997 for an alternative ap- 
proach). They describe the polarization field in 
terms of two quantities scalar under rotation, 
E(n) and B{h) (in analogy with electric and mag- 
netic fields' properties under parity transforma- 
tion). The relation between E, B and Q, U is 
non-local. In fact, to construct scalars under ro- 
tation at point h, we need to average the Stokes 
parameters around circles centered at n (Zaldar- 
riaga 2001). From the harmonic expansions of the 
electric and magnetic modes one defines the po- 
larization spectra Cei and Cbi- 

Since the analysis presented in this paper con- 
siders only small patches of the sky, we can locally 
approximate the sphere as a plane and, hence, re- 
place the spherical harmonic decomposition with a 
plane wave expansion (the standard Fourier trans- 
forms) (Seljak 1997). This allows us to consider- 
ably simplify the calculations. 

In the flat sky limit, we can fix a common co- 
ordinate system and use this to define the Stokes 
parameters at every point 9 in the plane of the sky. 
In this way, we can directly compute the Fourier 
components of Q and U, 



Q{i) = 



U(£) 



j &6e- llS PI{6) cos(20 e -) (2) 

fa 
j* 



'>/(«) sin(2fc) (3) 



where <f>g is the angle of which the common axis 
must be rotated in order to have U = in the point 
9. Using the Fourier transforms one can estimate 
the polarization power spectrum through: 



(4) 
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with X = Q,U. 

In the small scale limit the relation between 
electric and magnetic modes and the Q, U Stokes 
parameters assumes a very simple expression, con- 
sisting of a rotation in the ^-space: 

E(£) = Q{l)cos{2(j} i ) + U{i)sm{2(j) i ) 
B(£) = -Q{£) sin(2^-) + U{£) cos(20 f ) . (5) 

The Ce£ and Cbi spectra are then computed using 
eq. 4. 

The same expressions used to define the CMB 
polarization spectra can be safely applied to fore- 
ground components and, in particular, to the syn- 
chrotron emission. However, we would like to 
stress some differences. First of all, the spatial 
distribution of synchrotron radiation, both in the 
total intensity and in its polarized components, ex- 
hibits strong non-Gaussian features and may de- 
pend on the sky position. For instance, the Galac- 
tic plane, the North Galactic Spur and discrete 
sources generally present an emission level that is 
much higher than the other parts of the sky. Now, 
considering the Test Region, Fig. 1 shows, to- 
gether with an extended and bright region, areas 
with extremely low intensity levels. Such strong 
variations in the polarization intensity are due not 
only to the different intrinsic level of polarization, 
which we expect to be more uniform on the Galac- 
tic plane, but also reflect the different importance 
of the depolarizing mechanisms that occur inside 
the emitting regions or along the line of sight. 
Moreover, the polarized emission shows a weaker 
dependence with respect to b than the total inten- 
sity. So, the polarization APS computed from this 
survey does not describe only the intrinsic spatial 
distribution of polarized emission and can locally 
vary according to the area considered. However, 
when we consider patches of the sky that are large 
with respect to the dimensions of the observed 
structures, the angular spectrum is expected to 
adequately describe the average spatial properties 
of the emission. In our analysis, we consider also 
several small patches (1° x 1°) in order to study 
how the spectrum shape changes among regions 
with faint and bright polarized emission. 

Differences between CMB and foregrounds con- 
cern also the B-mode. For cosmological signals 
induced by scalar perturbations this quantity van- 
ishes. It arises only in presence of vector or ten- 



sor perturbations and in any case its intensity 
level is smaller than the E-mode one (Zaldar- 
riaga & Seljak 1997). On the contrary, the syn- 
chrotron emission contributes, on average, with 
the same amount to both quantities and their 
spectra present a similar shape (Seljak 1997). 

In order to estimate the APS from observational 
data, we need to consider in the above expression 
(4) the characteristics of the instrument. If X{9) is 
the quantity measured by ATCA experiment on a 
square grid of N pixels covering a solid angle Cl, in 
the small scale limit the estimator of the angular 
power spectrum can be computed as 

C ^ = w t ? I(/)r(/)6 / 2 w ~* ' (6) 

I 

where X{£) = N^ 1 e~'^'®X(9) is the discrete 
Fourier transform of the data and corresponds 
to the quantities Q, U, E, B. The sums arc 
performed over the Ng independent modes with 
wavevector magnitude around £. In the present 
case the term b^ takes into account that the 
ATCA images were subsequently smoothed with 
an elliptical Gaussian window function of FWHM 
a x 9b = 87" x 67". The expression in Fourier 
space is 

b?= expj-i/.M./j, (7) 

where the matrix M embodies the beam proper- 
ties. We checked that adopting a circular beam 
with 9 avc = 77" would make very little differ- 
ence. The contribution of the noise, which is 
smoothed by the window function b^, is subtracted 
by means of the factor w^ 1 = Slfccr 2 , i.e. the pixel- 
independent measure of noise (er is the rms noise 
amplitude, and fib — 8 a x 9b). 

Together with the electric and magnetic modes, 
in the following analysis we conside r the AP S also 
for the polarized intensity, PI = \J Q" 1 + U 2 . If Q 
and U are two Gaussian quantities, PI will have a 
Ricean distribution and its standard deviation will 
be api — 0.66 <tq,u- The polarized intensity is a 
scalar quantity and can be expanded in ordinary 
spherical harmonics; its power spectrum Cpn can 
be, therefore, defined in a way equivalent to the 
CMB temperature fluctuations spectrum. In the 
small scale limit, an expression equivalent to equa- 
tion (6) can be used to compute Cpn. In this case, 
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the contribution of noise (wp T ) will be lower than 
those of E and 5 by a factor ~ 0.43. The func- 
tion b£ is not exactly the antenna beam used for 
Q and U , however this can be considered a good 
approximation for b^. The differences with the ac- 
tual smoothing function become meaningful only 
at angular scales very close to the telescope reso- 
lution (they are responsible for the increasing of 
Cpn that we find at £ > 8000; see figures). 

The Cne spectrum must not be confused with 
the polarization power spectrum Cpg defined in 
Seljak (1997), where C P£ = C Qe + C ut = C E e + 
Cpe, and clearly has a different physical mean- 
ing with respect to Cpm. The latter, unlike the 
E- and B-modes, does not provide a complete 
description of the polarization field because it is 
related only to the intensity of the polarization 
without giving any information on its position an- 
gle. Only if the polarization angle is uniform in- 
side all the survey area, the equality Cpn = Cpi 
is ensured, while in general Cpn and Cpg are ex- 
pected to show different behaviors. This is surely 
the case for CMB polarization, as shown in Fig. 
2. There we report the -E-mode (solid line) and 
PI (dashed line) spectra resulting from Fourier 
analysis on CMB simulated maps. We chose a 
standard cosmological model with only scalar per- 
turbations, where therefore Cei = Cpe, and sim- 
ulation boxes of 10° x 10°. 

As discussed in the Section 1 , according to pre- 
vious analysis of the synchrotron maps by Brus- 
coli et al. (2002), the polarized intensity APS has 
a slightly faster decrease than the E and £?-modc 
APS at £ < 800. This is not surprising because we 
expect that the E and B modes vary more rapidly 
than PI, due to changes both in the intensity and 
in the angle of polarization. 

2.3. Results 

We study the APS of the synchrotron polariza- 
tion in several square patches of the survey with 
different dimensions. Because of the limited sky 
area considered, the Fourier approach is suitable; 
after computing the Fourier components of the 
data, Q{£), U(£) and PI(£), equation 6 is applied 
to estimate the spectra for the modes E, B and PI 
(the spectra for Q and U are redundant, as shown 
in the previous section). A standard least-squares 
method was used to fit each curve to a power law, 
Cxi = Ax(~ ax , in suitable ^-ranges (see Tucci et 



al. 2000). 

Fig. 1 shows the areas of the survey where the 
power spectra are computed: a 4° x 4° box cen- 
tered in (I, b) = (329°, 1?5), which is the largest 
square area that can be extracted from the survey; 
two 3° x 3° boxes covering nearly all the high- 
sensitivity part of survey; and six 1° x 1° boxes, 
in order to sample regions with different features 
in polarization. 

In the Test Region, 21 compact sources showing 
linear polarization were identified (see Table 1 in 
G01). Some of these can be readily seen in Fig. 
1. We tried to choose the 1° x 1° boxes in regions 
where the number of compact sources is not more 
than one. Also for the larger boxes we can fairly 
suppose that the number of compact sources is 
low enough not to affect the estimate of the power 
spectra. 

First of all, let us discuss the results for Cei and 
Cpe in the 4° x 4° area. In Fig. 3 we report only 
Cei (solid line); the results for B-mode are similar 
to those for i?-mode. In the range 600 < £ < 6000 
such spectra are well approximated by a power 
law with ole ^ ol b ~ 2.7 -r 2.9 (see Table 1). We 
compare the E-mode spectrum to that resulting 
from the observations in the same sky area by the 
Parkes telescope (D97). In the plot the Parkes 
spectrum has been scaled to the ATCA frequency, 
assuming the frequency spectrum of synchrotron 
emission T syn cx with spectral index (3 = —2.5 
and —3 (dashed lines). We compute the best fit 
for 100 < £ < 800 and find a power law index of 
ole = 1.78 ±0.18. This value is significantly lower 
than the best-fit a e found from the ATCA survey, 
i.e. 2.85 ±0.07 (dotted line). Moreover, the ampli- 
tudes of the Parkes and ATCA spectra show a dif- 
ference of nearly one order of magnitude at scales 
600 £ 1000. In this range, corresponding 
to angular dimensions between 20' and 10', both 
Parkes and ATCA are able to detect polarization 
signals. This gap, therefore, must be related to 
the different frequency at which the observations 
were made and cannot be attributed to our choice 
for the synchrotron spectral index. In fact, to re- 
move the discrepancy one requires a spectral index 
of— 5-; — 4, i.e., values too far from the average 
estimate between these frequencies ((3 — —2.8, see 
Platania et al. 1998), even in a very peculiar re- 
gion. 

Quite different results come from the analysis 
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of the polarized intensity spectrum. For the same 
4° x 4° survey patch (see Fig. 4) the Cpu is well 
approximated by a power law with api = 1.66 ± 
0.07 (the dotted line is the best-fit in 600 < I < 
6000). The PI spectrum shows a much flatter 
shape than Ce,bi, an d looks like the extension at 
high I of the curves obtained from 2.4 GHz Parkes 
data. 

Similar results are obtained from the analysis 
of smaller areas. In Fig. 5 we directly compare 
the electric and magnetic spectra with Cpu, com- 
puted on two 3° x 3° regions that sample all the 
low-noise part of the survey. Table 1 reports the 
best-fit parameters of the curves. These plots 
make evident how the E- and B-mode spectra 
have a similar shape to each other but separate 
from Cpu as they move toward low I. 

In Fig. 6 we plot Ce(. and Cpu obtained in 
1° x 1° regions, chosen between high- and low- 
polarized emission areas (see the rms values of the 
polarized intensity Firms in Table 2). We use let- 
ters a through / to label six areas. The best-fit 
parameters for the curves are reported in Table 
2. We found that the slope of both E-, i?-mode 
and PI spectra in 1° x 1° regions are quite inde- 
pendent of the area considered, with only small 
fluctuations around the values found in the bigger 
regions. For example, the spectral slope does not 
significantly differ between the regions a - b, that 
are within the brightest area of the Test Region, 
and e - /, where only faint polarized emission is 
observed. The only exception is the d case, which 
corresponds to a peculiar region inside which 1.4 
and 2.4 GHz data show anti-correlation, probably 
produced by internal Faraday depolarization oc- 
curring at 1.4 GHz. In this area Ce.bi an d Cpu 
are much steeper than in the other regions, indi- 
cating that internal depolarization erases intrinsic 
small-scale structures. 

2.4. ISM structures, Faraday screens and 
point sources 

The strong differences between Cpu and Ce,b£ 
found at 1.4 GHz but not at 2.4 GHz, as well as 
the amplitude gap in E-mode APS shown in Fig. 
3 can be interpreted if we assume that polariza- 
tion structures are not intrinsic to the emitting re- 
gions. As already discussed by G01 and in section 
2.1, an alternative mechanism to produce small- 
scale structures is the Faraday rotation along the 



line of sight due to foreground screens. In detail, 
if a smooth polarized background emission passes 
through an interstellar medium (ISM) with spa- 
tial variations in the rotation measure, the linearly 
polarized radiation is affected by a Faraday rota- 
tion of the polarization angle that varies spatially 
with the line of sight. Hence, if the angular di- 
mension of these ISM structures are sampled by 
the interferometer, the Faraday rotation will in- 
duce detectable variations in Q and U, and the 
E and B spectra will gain extra power on these 
scales. Visual inspection of Fig. 5 of G01, report- 
ing the polarization angle, apparently shows many 
sharp-edge structures with size of the order of 10', 
supporting our interpretation. 

The situation is quite different for the polar- 
ized intensity spectrum, Cpig. The PI value is 
not changed by a rotation of the polarization an- 
gle, but it could be affected by depolarization ef- 
fects that the rotation induces. Some areas, where 
Faraday rotation may be the cause of a strong re- 
duction of the polarized intensity, are found by 
G01: for instance, the existence of two large re- 
gions devoid of polarized emission is explained by 
beam depolarization due to large RM variations 
on small scales. Their position is, however, at the 
edges of the survey and their contribution to the 
APS of our selected regions can be considered neg- 
ligible. Internal Faraday depolarization could be 
the cause of the anti-correlation of the polarized 
emission between 1.4-GHz and 2.4-GHz data, ob- 
served in the "spur" region (labeled as d, cfr. Sec- 
tion 2.3). However, apart from this peculiar re- 
gion, the low RM values measured in the rest of 
the ATCA survey suggest that Faraday depolar- 
ization should not affect our estimates of the po- 
larized intensity spectrum, especially in the large 
boxes. This is confirmed by api values found in 
1° x 1° regions, which result weakly dependent on 
the area considered and on the magnitude of RM 
measured there (see Table 2). So, in the absence of 
depolarization, the Faraday rotation changes only 
the polarization angle without affecting the inten- 
sity and consequently Cpie- 

The effects of Faraday screens on the polar- 
ization emission have been discussed by previous 
authors. Wieringa et al. (1993) were the first 
to notice that the polarized component of the 
galactic background at 325 MHz is characterized 
by a patchy pattern on arcminute scales. They 
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detected features like narrow filaments of 5'— 10' 
width and "clouds" with angular sizes of some tens 
of arcminutes, across which there are only small 
changes in the polarization angle. These struc- 
tures were found at different latitudes and were 
just interpreted in terms of Faraday modulation 
by foreground ionized gas clumps. More recent 
interferometric observations (Duncan et al. 1998; 
Grayetal. 1999; Haverkornet al. 2000) confirmed 
these features in polarization even at 1.4 GHz. 

In Faraday-induced polarization structures ob- 
served by Wieringa et al. (1993), excess RMs from 
compact screens were measured with typical mag- 
nitude ~ 5 radm -2 , corresponding to A<f> ~ 13° 
at 1.4 GHz. Although small, these RM values suf- 
fice to produce structures which are detectable by 
the ATCA. Following the appendix in G01, if uni- 
formly polarized emission with intensity Pq passes 
through a compact cloud of rotation measure R c , 
the interferometer detects polarized emission in 
the direction of the cloud with apparent intensity 

P dct = 2P sin(.R c A 2 ), (8) 

and measures the RM value of the cloud as R c /2. 
Then, using R c = 10 radm -2 , at 1.4 GHz, the 
polarization induced by a compact cloud will be 
close to 100% of the background component. At 
2.4 GHz (the frequency of the D97 survey) this 
percentage is reduced by about a factor of three, 
without considering the depolarization effects due 
to the lower resolution. 

The above discussion assumes that none of the 
features of angular spectra should be attributed to 
extragalactic point sources. This is consistent with 
estimates: from Tegmark and Efstathiou (1996), 
using a VLA sample of 1.5 GHz sources, flat inten- 
sity spectra Cf/ = (3 13) x 10~ 8 K 2 are derived 
for limiting fluxes in the range (0.1 -f- 1) Jy. From 
Toffolatti et al. (1998) considering an evolution 
model for galaxies, for the same limiting fluxes 
and the (frequency) spectral index in the range 
-(2 -7- 2.3) we get Cf e s = (0.2 + 6.6) x lO" 8 K 2 . A 
more conservative upper limit can be derived from 
our own analysis of intensity spectra in U99 (Tucci 
et al. 2001, Bruscoli et al. 2002). Some patches 
exhibit very flat spectra there, a/ ~ 0; assuming 
that they are dominated from point sources, we 
get the limit Cff < 5 x 10~ 7 K 2 . From this num- 
ber, adopting a radio-source polarization degree of 
5% [in agreement with De Zotti et al. (2000)] we 



get C% s e < 1.3 x 10~ 9 K 2 for X = P, PI, and 
we conclude that the contribution of point sources 
should be negligible in the whole range £ < 6000. 
This is certainly consistent with the simple behav- 
ior of the angular spectra that we find. 

3. Discussion 

In this paper, for the first time, we extend the 
study of the angular power spectrum for the polar- 
ized component of the Galactic synchrotron emis- 
sion to arcminute scales, i.e. up to t ~ 10 4 . To 
reach such scales we needed high-resolution data, 
which were provided by the ATCA observations of 
a small patch of the Galactic Plane at 1.4 GHz. 

In the paper we compute the polarization spec- 
tra for "electric" and "magnetic" modes, plus the 
spectrum for the polarized intensity. We find that, 
in the range 600 < I < 6000, both C E t and 
C bi can be well approximated by power laws with 
slopes as — as ~ 2.7-r- 2.9. Such spectra are 
significantly steeper than those arising at I < 800 
from low-resolution data. Moreover, their ampli- 
tude, if compared to the spectra obtained by the 
Parkes telescope in the same sky area, turn out 
to be higher by nearly one order of magnitude at 
angular scales between 20' and 10'. These pecu- 
liar behaviors are well interpreted as due to the 
small-scale modulation of a relatively uniform po- 
larized background by Faraday rotation along the 
line of sight. On the contrary, we believe that our 
estimates of Cne, whose slope (c*pj ~ 1.7) is in 
agreement with D97 data at 2.4 GHz, are not af- 
fected by Faraday effects and fairly describe the 
intrinsic spatial distribution of the polarized emis- 
sion. 

An interesting point, which arises from our 
analysis, regards the distinctive meaning of Cpn 
with respect to Ce,bi- As we have discussed in 
section 2.2, PI is a scalar quantity and refers only 
to the intensity of the polarization without any 
information on its direction. We then expect that 
the APS for E- and -B-modes, that provide a com- 
plete description of the polarization field, do not 
have the same shape as the PI spectrum. The 
differences should be greater when the direction 
of polarization changes very rapidly. Deviations 
between Cei and Cpn are found in the CMB (see 
the results of the simulations in Fig. 2): these are 
not unexpected, because of the geometry of the po- 
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larization angle in the E-mode spots. In the case 
of synchrotron emission, the polarization direction 
for the diffuse component is quite smooth on large 
scales following the Galactic magnetic field. From 
low-resolution surveys the estimates of Ce,bi and 
Cne give only moderate differences in the spec- 
tral shape (see Bruscoli et al. 2002). However, 
when small scales are considered, fluctuations in 
magnetic fields, discrete sources and also Faraday 
effects contribute to amplify the variations in the 
polarization spectra, as the present results high- 
light. 

The extrapolation of our results for Cei and 
Cbi to higher frequencies should not be regarded 
as a trivial matter, since Faraday effects are sub- 
stantial at 1.4 GHz while they become negligible 
at a few tens of GHz. The electric and magnetic 
spectra are strongly affected by the Faraday rota- 
tion along the line of sight, showing a steep slope 
(ctE.B ~ 2.8). The polarized intensity spectrum, 
instead, can be more reliably extrapolated to the 
"cosmological" frequencies, because, as discussed 
in section 2.4, it is not affected by Faraday rota- 
tion, except in severe cases when significant de- 
polarization is occurring. The power index api 
is less than 2 independent of the region analysed, 
with a value of 1.66 ± 0.05 in the 4° x 4° box. 

We have seen that analyses on the synchrotron 
polarization spectrum in the literature indicate a 
moderate slope (ax ^> 2 with X = E,B, PI) on 
angular scales I < I0 3 . Now ruling out a signifi- 
cant contribution from point sources, we confirm 
this result also at t ^ 10 4 for Cpn, and we show 
that synchrotron emission is rather rich in small 
scale structures. Hence, contrary to what is usu- 
ally assumed, it might be a relevant contaminant 
in CMB polarization measurements at very small 
scales. It remains to be seen if the APS as deduced 
by ATCA observations of the Galactic plane is a 
common feature in regions at high galactic lati- 
tudes. In general we expect the polarized syn- 
chrotron emission to be fainter in regions far out 
of the Galactic plane, except in very bright ar- 
eas. In one such region Bruscoli et al. (2002) 
estimate Cxe for I < 100, finding a spectrum be- 
havior consistent with those of the Galactic plane 
(1 < ax < 2). Even if very bright regions are 
not typical at high latitudes, these play an impor- 
tant role in the process of foreground subtraction 
in CMB experiments; this is why we have to put 



great care into the study of the APS of these re- 
gions. 
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Fig. 1. — Linearly polarized intensity, PI = 
\/Q 2 + U 2 , in the Test Region. The image is ob- 
tained combining all nine spectral channels. The 
intensity scale runs from 0.4 to 9.5 mJy beam . 
The square boxes indicate the areas where we com- 
puted the angular power spectra. 



This 2-column preprint was prepared with the AAS IATfrjX 
macros v5.0. 
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Table 1: Best -fit parameters for angular power spectra 



Survey 1 ' 


Box 


Centre 


A E 


a E 


a B 


An 


Oipi 


£-range 




(deg) 


(l,b) 


(K 2 ) 






(K 2 ) 




(xlO 3 ) 


D97 


5x5 


(329,1.5) 


9.1xl0~ 4 


1.78 ±0.18 


1.62 ±0.19 


3.6xl0~ 4 


1.68 ±0.30 


(0.1,0.8) 


G01 


4x4 


(329,1.5) 


93.1 


2.85 ±0.07 


2.74 ±0.06 


4.2xl0 -3 


1.66 ±0.07 


(0.6,6) 




3x3 


(330.5,1.5) 


60.2 


2.76 ±0.09 


2.65 ±0.10 


6.0xl0~ 3 


1.67 ±0.08 


(0.6,6) 




3x3 


(327.5,1.5) 


63.0 


2.84 ±0.09 


2.70 ±0.10 


3.7xl0" 3 


1.68 ±0.08 


(0.6,6) 



1 For the D97 survey the A x values are at 2.4 GHz, for the Test Region (G01) at 1.4 GHz. In the last 
column we report the range of angular scales considered in the fit. 



Table 2: The best-fit parameters for 1° x 1° areas in the €-rangc [1000, 6000] 



Box Label 1 


Box Centre 


PT 

± J rms 


A E 


a E 


as 


Api 




api 




(/, b) 


(K) 


(K 2 ) 






(K 2 ) 






a 


(330.6,1.0) 


0.46 


168 


2.84 ±0.25 


2.58 ±0.24 


2.1x10" 


-2 


1.79 ±0.20 


b 


(329.7,1.1) 


0.54 


2.9xl0 3 


3.18±0.26 


3.24 ±0.30 


0.89 




2.24 ±0.22 


c 


(327.7,1.1) 


0.38 


137 


2.91 ±0.30 


2.97 ±0.28 


3.4x10" 


-2 


1.93 ±0.27 


d 


(328.5,2.6) 


0.38 


1.8xl0 5 


3.84 ±0.34 


3.53 ±0.35 


5.1 




2.58 ±0.29 


e 


(327.0,2.3) 


0.26 


616 


3.18±0.35 


3.07 ±0.32 


1.0x10" 


-2 


1.88 ±0.38 


f 


(330.4,2.3) 


0.31 


233 


3.01 ±0.38 


2.64 ±0.39 


4.4x10" 


-3 


1.68 ±0.35 



t The letters correspond to the panels in Fig. 6. 
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Fig. 2. — The CMB spectra for E-mode (dotted line is the input spectrum, the solid line is the computed 
spectrum from Fourier analysis) and for the PI field (dashed line), from a CDM model with a secondary 
ionization optical depth r, „ = 0.2. 
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Fig. 3. — E-mode spectrum from the 4° x 4° region (solid lines), compared to the result from D97 data 
(dashed lines) in a 5° x 5° area centered on the same position. The dotted line is the best-fit power law in 
the ^-range between 600 and 6000 (see Table 1). 
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Fig. 4. — As in Fig. 3, but for polarization intensity spectra. 



15 



8 (arcmin) 
100 10 



1 



6 (arcmin) 
100 10 1 




Fig. 5. — Ce£ (solid lines), Cbi (dotted lines) and Cm (dashed lines) in two 3° x 3° boxes, whose centers 
are shown in the panels. 
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Fig. 6. — Ce£ (solid lines) and Cpn (dashed lines) in six 1° x 1° boxes. The labels in the plots identify the 
regions where the spectra are computed. 
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